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Spectroscopic Characterisation of 250/im-Selected 
Hyper-Luminous Star Forming Galaxies 
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ABSTRACT 

We present near-infrared spectroscopic observations from VLT ISAAC of thirteen 
250/im-luminous galaxies in the CDF-S, seven of which have confirmed redshifts which 
average to (z) = 2.0±0.4. Another two sources of the 13 have tentative z > 1 identi- 
fications. Eight of the nine redshifts were identified with Ha detection in H- and 
K-bands, three of which are confirmed redshifts from previous spectroscopic surveys. 
We use their near-IR spectra to measure Ha line widths and luminosities, which aver- 
age to 415±20kms- 1 and 3xl0 35 W (implying SFR ffQ -200M o yr- 1 ), both similar 
to the Ha properties of SMGs. Just like SMGs, 250/im-luminous galaxies have large 
Ha to far-infrared (FIR) extinction factors such that the Ha SFRs underestimate 
the FIR SFRs by ^8-80 times. Far-infrared photometric points from observed 24/im 
through 870/im are used to constrain the spectral energy distributions (SEDs) even 
though uncertainty caused by FIR confusion in the BLAST bands is significant. The 
population has a mean dust temperature of Td = 52±6K, emissivity /3 = 1.73±0.13, 
and FIR luminosity Lfir = 3x10 13 Lq. Although selection at 250/mi allows for the 
detection of much hotter dust dominated HyLIRGs than SMG selection (at 850/xm), 
we do not find any £ 60 K 'hot-dust' HyLIRGs. We have shown that near-infrared 
spectroscopy combined with good photometric redshifts is an efficient way to spectro- 
scopically identify and characterise these rare, extreme systems, hundreds of which are 
being discovered by the newest generation of IR observatories including the Herschel 
Space Observatory. 
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1 INTRODUCTION 

Submillimetre Galaxies (SMGs) contribute significantly to 
the rapid buildup of stellar mass in the Universe at z ~ 2. 
However, their selection at 850/mi is inherent l y biased to- 



ward s colder-dust sourc es |E alcs ct al 
20041 ). Recent work (e.g. IChapman et alj 



2000; Blain ct al 



2004 ICasev et al 



2009c) has demonstrated that 850/im-faint, high-redshift 
Ultraluminous Infrared Galaxies (ULIRGs) exist and may 
contribute significant ly to the co s mic st ar formation rate 
density at its peak. I Casey et al.l |2009cf ) describe a pop- 
ulation of 70/im luminous galaxies at z ~1.5 whose in- 
frared luminosities exceed ~10 12 L© but are 850/im-faint 



Hubble Fellow; cmcasey@ifa.hawaii.edu 



due to hotter characteristic dust temperatures. Sparse in- 
frared data, particularly in the 50-500/tm wavelength range, 
along with poor volume density constraints have limited the 
interpretation of these submm-faint ULIRGs. Similar studies 
of other infrared-luminous galaxy populations, selected at 
24/im, 350/im or 1.2mm for example, present even more ev- 
idence for diverse populations of luminous, dusty starbursts 
at z Si 1 which do not necessarily intersect (see iDev et all 
120081 ; iBussmann et al.ll2009l ; lYounger et alj|2009h . 

The arrival of new infrared (IR) instruments— including 
BLAST ( the Balloon-borne Large-Aperture Submillimeter 
Telescope IPascale et al.ll2008r i, SCUBA2, LABOCA, and the 
Herschel Space Observatory— has opened up more exten- 
sive studies of these distant starbursts. BLAST'S deep map- 
ping of the Extended Chandra Deep Field South (ECDF- 
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Figure 1. The rcdshift distribution of BLAST 250/mi sources, 
including the low-redshift spectroscopically identified sources in 
D10 and 110 (gray points/histogram). Our spectroscopic sam- 
ple has a mean redshift 2.0±0.4 (black crosses, lined histogram), 
which fills in much of the high-rcdshift subsample of 250/im-bright 
objects. The distribution in photometric redshifts is illustrated as 
the dashed line histogram and are good to ~12%. The light gray 
vertical stripe blocks out the redshift regime where near-IR spec- 
troscopy cannot confirm redshifts since Ho falls between H and 
K bands. 

S) at 250/im, 350/im, and 500/im has, for the first time, 
led to z^2 ULIRG sele ction near the peak of their SED. 
Ilvison et all (|2010l l and iDunlop et ail |2009l ) describe the 
selection of these 250/im sources, along with their radio and 
24/im counterparts, in detail and match sources to photo- 
metric redshifts derived from the extensive ECDF-S multi- 
wavelength data. We also make use of longer wavelength con- 
straints from the LABOCA 870/im survey of the ECDF-S 
jWeifi et al.ll2009l V While most low redshift (z < 0.8) 250/im 
sources have spectroscopic identifications, none of the sus- 
pected high-redshift sources had spectroscopic redshifts. 

This paper presents new VLT ISAAC spectro- 
scopic observations of thirteen BLAST 250/im sources 
with Zphot > 1. With spectroscopic redshifts, we constrain 
the FIR dust SED (implying that they are HyLIRGs 
with Lfik ^ 10 13 Lq), measure dust temperatures, black- 
body emissivity, FIR luminosities, Ha luminosities and 
AGN/metal lines, and constrain the FIR/radio correlation 
for these high-z galaxies. Thr oughout we use a, Acdm cos- 
mology (jHinshaw et alJ [2009) with H =7l kms~ 1 Mpc" 1 
and n = 0.27. 



2 OBSERVATIONS & RESULTS 

Longslit spectroscopic observations of 250/im sources were 
obtained in November 2009 on the Very Large Telescope 
(VLT) Infrared Spectrometer And Array Camera (ISAAC) 
under excellent seeing conditions (0.3-0.7" in K- band). Spec- 
trosco pic c andidates wer e chos en from the Ilvison et al.l 
j20ld ) and IDunlop et all (|2009r ) (hereafter 110 and D10) 
ECDF-S BLAST Deep map samples (RMS sensitivity 



0250 = 11 mjy) with photometric redshifts above z —1 or 
undefined photometric redshifts (the latter caused by a lack 
of high quality photometry). 110 selected sources at >5<r, 
having folded in the confusion noise (~21mJy), resulting in 
flux densities S250 >59mJy. D10 selected sources at >3u 
without accounting for the confusion noise, so their source 
list has 5250 >33mJy. All had reliable 24/im and/or radio 
counterparts, which were then matched to K-band sources 
(with offsets < 1" to the radio/24/i m centroid) in archival 
MUSYC data jGawiser et alj|2006h for VLT spectroscopic 
targeting. The resulting candidate object list contained 20 
sources. 

We refer the reader to D10 and 110 for the analysis of 
the BLAST 250/im CDF-S map and source selection, as well 
as some source properties derived from ancillary data. The 
two papers present different detection thresholds (which are 
discussed in more detail as they relate to source density es- 
timations in section [3~Tj) and also use different counterpart 
identification methods to identify sources for photometric 
redshift fitting. D10 uses 24/im counterpart matching, while 
110 uses 1.4 GHz radio matching. In general radio match- 
ing is muc h more reliable ass uming the FIR/radio correla- 
tion holds (|Helou et al.| [l985) as there are far fewer poten- 
tial counterparts and good reason to suspect a FIR-bright 
source is also radio-bright. Identification at 24/im is a less re- 
liable alternative due to the large IR beamsize and density 
of sources. We caution the reader that our K-band iden- 
tifications are nearest neighbours to the radio and 24/im 
counterpart astrometry of D10 and 110 and that there is 
minor potential for misidentified counterparts. The offsets 
between 250/im peaks and K-band sources (which are effec- 
tively equivalent to the radio/24/im positions) range from 1- 
16", averaging to about 7", which is well within the beamsize 
of 250/im observations, however it is not possible to know 
if the counterparts have been identif ied correctly withou t 
high-resolution FIR observations (e.g. lYounger et "al"1l2010l ). 

We observed 13 of the 20 z p hot > 1 250/im sources 
searching for Ha or [OIII] in J, H, and K bands. The band 
of observations was primarily chosen based on the galaxies' 
photometric redshifts (where galaxies with z p i lo t >2 were 
observed in K-band and at z p hot <1.8 in H-band, and in 
J-band for the intermediate region). Galaxies were observed 
individually under varying seeing conditions which ranged 
from 0.3-0.7" seeing in K-band. Since all galaxies here are 
assumed to be unresolved, we varied the slit width according 
to seeing conditions, minimising it when possible to reduce 
sky line contamination. Galaxies were centred on the 2' long 
slit and observed in ABBA nodding mode with a 15" nod. 
On occasion, two candidates were within a 2' separation and 
placed on the same slit, with the maximum possible nod 
distance, which was sometimes 5-10". Data reduction was 
completed with ESO software combined with IRAF and our 
own IDL-based routines to obtain ID and 2D wavelength 
calibrated spectra. 

Only 13 of the sources were observed due to telescope 
time constraints. Nine of the 13 have spectroscopic red- 
shifts, seven of which are secure. Three of these seven sources 
(J033246, J033152, and J033243) were already identified in 
previous spectroscopic surveys using the Very Large Tele- 
scope/FOcal Reducer and low dispersion Spectrograph 2 
and the Gemini Near-In frared Spectrograph |Vanzella et al.l 
120081 ; iKriek et all l2008h at z =1.382, 2.336 and 2.122 re- 
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spectively. Our measured redshifts confirm these observa- 
tions. Six of the seven secure redshifts and the two ten- 
tative redshifts were measured from Ha detection (at >4u 
Ha signal-to- noise) . One of the eight Ha redshifts (J033129) 
would nominally be tentative, but it was spectroscopically 
identified in the rest-UV independently; the ninth secure 
redshift which is not based on Ha (J033151) has absorption 
features in the K-band at ~1.599 which agree with a rest-UV 
redshift of 1.605 obtained independently (Swinbank, private 
communication). The lines were identified as Ha using a 
combination of photometric redshift consistency and a lack 
of other line features (which would instead identify the line 
as [Oil] or [OIII], in the case of [Nil] or SII detection). 
The four sources which were not identifiable in emission ei- 
ther have very weak emission features or lie at redshifts in 
the range 1.7 < z < 2.0; Ha at these redshifts falls between 
H- and K- bands and is thus not detectable with near-IR 
spectroscopy. 

The nine galaxies have a mean redshift z = 2.0±0.4, 
and their redshift distribution, with respect to other 250/im 
sources, is shown in Figure [T] We determine that their pho- 
tometric redshifts (derived in D10 and in Rafferty etal., in 
prep, for 110 sources) are good to dz/(l + z) 12%, and 
we emphasise that this relatively small error implies that 
near-IR spectroscopic followup for ULIRGs with good pho- 
tometric redshifts is efficient. Table [T] summarises the galax- 
ies properties and Figure [2] shows their ISAAC spectra for 
those which were spectroscopically identified. The galaxies' 
names are derived from their positions in the K band. 

The spectra in Figure [2] are framed around the Ha 
emission for every source except J033151. Regions where 
emission lines in the sky's infrared spectrunQ are significant 
(with flux densities in excess of ~57s _1 nm -1 arcsec -2 m -2 
where 7 represents photons) are masked out in both ID and 
2D spectral renditions. The width of these sky lines varied 
according to the slit width of each observation, which varied 
from 0.3-0.8". 

We use the publicly released s moothed maps of the 
BLAST ECDF-S jPevlin et al.ll2009T ) to derive 350/im and 
500/iin flux densities for the D10 sample (both 350/im and 
500/im flux densities are published for the 110 sample). We 
use the K-band astrometry (where the VLT slit was placed) 
as position priors to extract BLAST flux densities since the 
maps are dominated by confusion noise. We note however, 
that because the maps are highly confused, they are likely 
to suffer from flux boosting on in dividual sources, espe cially 
for the fainter 250/im sources; see lChapin et all ((2010) for a 
detailed analysis of source extraction, confusion limitations 
with BLAST data, and SED fitting. 

The fluxes in Table [1] are not corrected for deboosting, 
but we do apply corrections when fitting FIR SEDs and de- 
riving Lfir, Tdust and j3. We use the signal-to-noise ratio 
of the 250/im detections to determine the deboosting fac- 
tor which corrects for the Eddingt o n bia s and boosting by 
confusion noise, as in lEales et al.1 |2009l ), figure A2. This 
results in flux densities which are ~55-75% of the original 
measured values for the BLAST bands given in Table [T] 



1 See the Gemini Observatories IR Background Spectra page 
(http://www.gemini.edu/?q=node/10787 ) for example sky spec- 
tra. 



We do not correct 870/im LABOCA points for deboosting 
since the 870/im flux densities are less likely to be boosted 
by flux from adjacent sources. This is because the sources 
which boost a high-redshift galaxy's 250-500/im flux density 
are far more likely to sit at lower redshifts than at higher 
redshifts (where the surface density of sources is more rare). 
These low redshift sources are also unlikely to be bright at 
870/im given typical blackbody SED shapes at low-2. We 
recognise that our application of a deboosting factor to the 
BLAST bands and not the LABOCA data is a simplification 
of a complex issue, but we consider it the most realistic con- 
straint on the FIR flux densities given the data which are 
available. We do not deboost the Spitzer MIPS photometry 
at shorter wavelengths since the source surface density at 
~70/im is similar to the source density at £ 200/im, yet the 
relative beam size is much smaller at shorter wavelengths. 

We choose not to propagate the uncertainty in the 
FIR deboosting factor (which is estimated to be as large 
as 50% at low 250/im S/N) into the SED uncertainty be- 
cause the deboosting factor is not independent between 
bands. Despite the large uncertainty, the FIR flux den- 
sities for individual sources would be deboosted by simi- 
lar factors. A correlated deboosting factor between bands 
would imply less uncertainty in derived Td us t or /3 than 
blindly adopting the deboost uncertainties. We test the cor- 
relation of the deboosting factor by using the FIR colours 
(Saso/ff aso and S250/S500) o f Herschel SPIRE and PACS 
sources (|Amblard et alj|20ich . Using Monte Carlo tests, we 
remove the contribution of a single potential boosting source 
by subtraction of an arbitrary 250/im flux and the associ- 
ated 350/im and 500/im fluxes associa ted with the colours o f 
a randomly selected galaxy from the lAmblard et all (|2010t ) 
sample. We find that the FIR luminosity does not vary by 
more than ±0.1 dex and that dust temperature varies by 
about ±9K. If two contaminating boosting sources are in- 
corporated with different FIR colour properties, the variance 
on the fitted Lfir and Tdust decreases further. We discuss 
the impact that the deboosting factor uncertainty has on 
our final conclusions more in § 13.51 

We measure 870/im flux densities (at the K-band po- 
sition s) from the LABOCA map of ECDF-S jWeifi et all 
2009h. Six source s have 5870 ~ 4mJy and are listed in 
Weifi et all (|2009l ). J033246 is claimed as an SMC in D10, 
but its 870/im-peak is >30" away from its K-band position. 
This implies that 7/13 (~53%) of our sample are submm- 
faint and would be excluded from traditional SMC surveys. 
All galaxies except J033212 are also radio detected in VLA 
data at >30 /iJy. 

The galaxies' rest-frame near-IR photometry is con- 
sistent with stellar emission, from which we derive stellar 
masses from rest- frame J/-band (Table Q , using the meth- 
ods described by lHainline et al. (2009). Measuring the ab- 
solute magnitude of a galaxy near its 1.6/im "stellar bump" 
provides the most accurate measure of its stellar mass, how- 
ever it is reliant on the assumption of a constant mass to 
light ratio (here we assume M/L = 3.2), reddening proper- 
ties, and minimal AGN contribution to near-IR flux. The 
stellar mass estimates are uncertain by ~0.3dex. The near- 
IR photometry is also be used to infer AGN content, since 
a flux excess at 8/im (significantly above stellar population 
model fits) is indicative of power law emission from an AGN. 
None of our sources have >2cr 8/im flux excesses. 
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Figure 2. The VLT ISAAC spectroscopy of nine 250/im BLAST sources around their Ha emission (0.63-0.69/xm rest wavelength), 
except in the case of J033151 when no Ha observations were obtained. The first seven sources have secure redshifts, the next two are 
tentative, and the last is a stacked Ha spectrum for the six secure Ha sources. Observations were tak^r^jjjjjHpgr^jK j!^^|^v\^$(-j^^puf*h 
and Id projections are shown for clarity; scaling is optimised for viewing spectral features which are marked by vertical dashed lincsTe.g. 
Ha, Nil). The 2d spectra show an angular scale of 12" top to bottom and the Id spectra are extracted within 0.6". Sky line emission 
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Table X. Multi- wavelength properties of BLAST 250 /^m Galaxies 
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ID 6 


2spcc 


^phot 


S'24 
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S 3 so d 


S 5 oo d 


•Si. 4 


5 8 70 d 
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P QlR f 










OJy) 


(mjy) 


(mjy) 


(mjy) 


0*Jy) 


(mjy) 




(10 13 L ) 


(K) 
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J033129. 874-275722. 40 


J033129 


1.482 


1.57 


270 


91.6±11.0 


54.4±8.7 


46.3±6.2 


144±16 


5. Oil. 5 


SMG 




45.9±3.9 


1.2±0.3 2.7±0.5 


J033151. 088-274436.91 


J033151 


1.599 


1.91 


520 


74.0±10.8 
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2.342 c 


2.30 
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2.09 
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30.1±10.9 


32.4±8.5 
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5.7±1.4 


SMG 
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0.8±0.7 1.6±0.4 


Trn^94fi 39Q 97^97 01 
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1.37 
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9c -1 4. inn 
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-1.3il.4 


SFRG 
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33.2±8.5 
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Table Notes. a Galaxies are split into three categories: 'detections' (sources which have reliable redshift identifications), 'tentative' 
(poor quality redshift identifications), or 'non-detections' (no visible emission features). All redshift identifications are based on Ha 
detection except J033151. The galaxies with 'tentative' identifications are included in all figures and tables of this paper but are 
excluded from the primary analysis points in section [3] so as not to affect the interpretation of this paper. 
b ID is the identification of the 250/^m source taken from D10 or 110. Those from 110 are of the form J033XXX and correspond to the 
first half of the BLAST name given in table 1 of 110. Those from D10 appear as two to four digit numbers and can be found as the 

BLAST IDs in table 1 of D10. 
c Three sources have confirmed redshifts from previous spectroscopic surveys llVanzella et alJl200St iKriek et al.1l200Sl) . 
d S350, S50O1 and 5870 are measured directly from BLAST/LABOCA ECDF-S maps at their respective wavelengths, using the K-band 
astrometry for sources in the D10 sample, and 24/jm and 1.4 GHz flux densities are based on nearest neighbour matching (described in 
110). The flux densities at 250/^m, 350/^m and 500/^m have not been corrected for flux boosting, and their uncertainties here only 
represent instrumental uncertainty; they should be combined in quadrature with the confusion noise (~21 mjy) for an accurate 

representation of flux uncertainty. 

e A galaxy's class is either SMG (su bmm galaxy ) or S FRG (submm-faint radio galaxy) based on its inclusion as a significant detection 

in the Weifi et a lj (1200^) sample, i.e. if its 870/^m flux density is ^ 4 mjy. 
* The ratio of IR luminosity to radio luminosity (as calculated in 110 using a = 0.75, see their section 2.2), qm, is given in the last 

column. 



Two of the 13 observed sources, J033151 and J033152, 
are X-ray detected above the luminosities which would cor- 
respond with their star formation rates (Lx ~ 10 44 ergs _1 ). 
Only one individual source shows obvious signs of containing 
a luminous AGN, from its radio flux excess and detection 
in the X-rays: J033152. From this AGN estimator (and the 
analysis in the ensuing section about Ha properties), we in- 
fer roughly that 20±15% of 250/im-bright sources have signs 
of dominant AGN. 

2.1 Ha Properties 

We measure Ha line widths and [Nil]/ Ha ratios in order 
to infer AGN content from the six BLAST sources for which 
we have secure Ha observations (i.e. not including J033151). 
After deconvolving the measured full width half maxima 
(FWHM) with the instrumental resolution measured from 
skylines in the vicinity of Ha (~6.5A in K-band and ~4.4A 
in H-band) we find that our Ha lines have an average rest- 
frame FWHM of 415±20 kms -1 and span the range 150- 
800 kms -1 (for the six galaxies with secure Ha detections). 
The individual Ha properties of each galaxy are given in 
Table [2] We plot the Ha FWHM against Ha luminosity 



in Figure [3] The Ha-inferred star formation rates average 
to 210±160 Mg yr" 1 , which requires a substantial extinction 
factor to account for the star formation observed in the FIR 
(on the order of 2000 M yr" 1 ). The mean SFRfir/SFR^ 
ratio for the sample is 36±22, w hich is comparable to the 
ratio for the SMG population of ISwinbank et all (|2004l ) of 
31±15. The subset of our sample which is submm-faint 
(SFRGs) also have similar SFR ratios, averaging 37±24. 

For the six galaxies which have secure Ha observa- 
tions, a stacked spectrum is shown in Figure [2] which we 
use to measure the aggregate line emission properties of 
the sample. The Ha line width of the stacked spectrum is 
530±280 kms" 1 (statistically indistinguishable from the in- 
dividual Ha measurements or the mean SMG line width, 
390 kms -1 ), and a line luminosity corresponding to a star 
formation rate of 190M©yr _1 . The line width is slightly 
larger than the mean line width for the sample likely due 
to signal-to-noise limitations of the original data. Both line 
width measurements, 415 kms -1 and 530 kms" 1 , are con- 
sistent with the dynamics of active star forming HII regions, 
except the high FWHM outlier: J033243 at z = 2.123 with 
FWHM = 800 kms" 1 . J033243 is also the second brightest 
Ha emitter with a high Ha implied SFR, 335M0yr _1 ; its 
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have Td and Lfir as free parameters. The advantage of al- 
lowing emissivity to vary in the first model allows a reassess- 
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Figure 3. The relation between Ha FWHM and line lumi- 
nosity is shown below , for the BLAST sample, SMGs (crosses; 
ISwinbank et al.1l2004l), an d normal z~2 galaxies (e.g. LBGs, gray 
triangles; lErb et al.ll2003h . The Ha line luminosity is converted 
to a SFR at right and the Ha FWHM is converted to a virial 
mass estimate (using the assumption of rw2 = 3 kpc) at top. The 
BLAST 250/tm sample has similarly broad and strong Ha emis- 
sion to the SMG population, and is brighter than the normal z~2 
galaxies. The two galaxies with tentative redshifts are shown as 
empty circles. This plot suggests that 250/jm-luminous galaxies 
might share intrinsic, physical properties with the SMG popula- 
tion. 



SFRFiR/SFRfla = 7, which is the lowest SFR ratio of the 
sample indicative of a less Ha obscuration. 

We convert the [Nil]/ H a ratios to 12 + l og(Q /H) us- 
ing the methods described bv lMaiolino et all (|2008h . How- 
ever two sources (J033249 and J033212) have O/H limits 
>9.25, which corresponds to very strong [Nil]/ Ha( ^ 0.5). 
The O/H and [Nil] / Ha indicators satur ate at metallici- 
ties above solar (see iPettini fc Pagelll2004l y and additional 
contribution from either AGN or shocked gas can in crease 
the [Nil]/ Ha ratio further (e.g. Ivan Dokkuml 12005V The 
remaining five secure detections have 12 + log(0/H) val- 
ues which average to 8.64T0.08. The measured [Nil]/ Ha 
ratio for the stacked Ha spectrum implies a metallicity of 
12 + log(0/H) = 8.75to ?o (in agreement with the average 
for the individual measurements). 



2.2 Dust SED Fitting and FIR/radio Correlation 

We fit the MIPS (70/^m, and 160/im, where available), 
BLAST (250/xm, 350A«n, and 500/im), and LABOCA 
(870/^m) flux densities to two different FIR dust models. For 
FIR SED fitting, we correct the BLAST flux densities for 
boosting by confusion noise as mentioned in the beginning 
of this section. Both FIR SED models assume a modified 
blackbody emission curve with a single dust temperature: 

S " ^ exp(hv/kT dust ) - 1 (1) 

where S v , the flux density, is a function of rest frequency v, 
the emissivity /3, dust temperature Td, and FIR luminosity 
Lfir (which governs the normalisation of the function). The 
first model allows /3 to vary (the "beta-free" model) while 
the second model fixes emissivity to /3 = 2. Both models 



ULIRGs in past studies (e .g. /3 



P =1.5 or 2.0; 


Chapman et al.l 


1 Younger et al.l 


2009). In addi- 



tion, our measurements of j3 are made independent of any 
a priori constraint on Td or Lfir- We choose to make the 
second model rigid as fits from the first model can be un- 
physical, as might be the case if the FIR flux densities are 
particularly faint or affected by source confusion. 

Only J033212 and J033237 are poorly fit to a beta-free 
model (these are the two galaxies with tentative redshift 
identifications), since they do not have 70/im data and have 
unconstraining upper limits in the FIR. We use only the 
fixed /3 = 2 model for these two. The remaining seven galax- 
ies have reliable beta-free SED fits, and from them we mea- 
sure /3, Tdust, and L F ir( 8 -looo^m) (summarised in TableQ]). 
Both fixed beta and beta- free fits are shown in Figure [4] We 
find a mean emissivity of /3 = 1.73±0.13 and a mean dust 
temperature of Td = 52±6 K. 

The FIR luminosities (8-1000/im) must be corrected to 
account for mid-infrared (8-25^m) emission from PAH and 
power law sources (e.g. iMenendez-Delmestre et all l2009h 
above the single FIR modified blackbody. We tether the 
iPope et alj i|200Sl) SMG SED to 24/xm flux densities (as seen 
in Figure |4j to estimate the luminosity deficit of the single 
temperature blackbody. This deficit varies substantially ob- 
ject to object due to the large spread in 24^im flux densities 
and blackbody properties in the 8-25^im Wein tail. On av- 
erage, we find that the contribution of the PAH and AGN 
emission account for 0.04T0.03 dex of luminosity which we 
add to the FIR luminosities as a correction factor. Although 
the mid-IR properties of the sample can vary substantially, 
this deficit translates to no more than a ~10% increase in 
FIR luminosity for these >10 13 Lq systems. The corrected 
luminosities are given in Table [U 

We also over plot the composite SMG spectrum, from 
I Pope et al.l (|2008T ), normalised to the integrated 24/im flux 
density in Figure [4] While the SMG composite is carefully 
derived based on mid-IR to FIR data of SMGs to date, it 
fails to fit the BLAST FIR data on a case by case basis. 
In some cases, it under/overestimates the FIR luminosities 
by Tldex. This illustrates how a 24/^m-normalised SED 
fitting procedure , which is common in the literature (e.g. 
iDesai et al.ll2009T ) places poor constraints on the breadth of 
FIR properties of ULIRG samples, especially in the absence 
of direct FIR measurements. Recen t high-resolution FIR ob- 
servations (e.g. lYounger et ail 2010) have demonstrated that 
24/im counterparts are often misidentifications and do not 
correspond to the FIR luminous source. 

Although multiple dust temperature blackbodies are 
found t o fit well to local ULIR Gs in the literature (for exam- 
ple, see lClements et al] [2010V strong assumptions must be 
made regarding the FIR luminosity or normalisation, to de- 
compose the sparse FIR data down into multiple blackbody 
components. Given the uncertainty of the FIR luminosities 
or flux densities at any given wavelength, we decide to forgo 
multiple dust temperature fitting for well constrained, sin- 
gle dust temperature blackbody fits. If multiple blackbodies 
provide a more physical SED fit, then our derived emissivi- 
ties, from the single blackbody fits, could be underestimated. 
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Table 2. Ha and stellar properties of BLAST 250 fim Galaxies 



NAME 


2 






FWHM Ha 


<[NII]/Ha) a (O/H)" SFR Ha 


SFR F i R 


SFR ra( jio 








(Wm" 2 ) 


(10 35 W) 


(kms" 1 ) 


(Mq yr" 1 ) 


(Moyr" 1 ) 


(Moyr" 1 ) 


(M ) 



J033129 


1.482 


4.1x10- 


19 


(s.6ir 8 ) 


190±60 


0.13 


8.60 


^±1 


21001"°° 


1300+||2 


lxlO 10 


J033151 


1.599 






(53+J) 








419+?° 


3800« 


1000+ 430 
±uuu_ 300 


3xl0 10 


J033152 


2.342 


1.3x10" 


18 


490±40 


0.74 


>9.25 


14000+™°° 


280001^°° 


3xl0 10 


J033204 


2.252 


5.2x10" 


19 


(20l 3 3 ) 


360±40 


0.14 


8.64 




6800!^ 


3300+ 1 """ 


5xl0 10 


J033243 


2.123 


1.3x10" 


18 


(42±|) 


800±50 


0.18 


8.71 


335 


2400lt 3 °° 


20001™° 


4xl0 10 


J033246 


1.382 


4.2x10" 


19 


(4-9tS:S) 


150±70 


0.13 


8.59 


3911 


6801?°° 


650+ 180 


lxlO 10 


J033249 


2.326 


5.8x10" 


19 


(25l 3 3 ) 


360±50 


0.14 


8.64 


194+*> 


1400011°°° 


6100+ 1 ™" 


3xl0 10 



J033212 


1.93 


3.5xl0" 19 


(9.3l^) 


310±50 


0.55 


>9.25 


7 q+i5 

' a -13 


3800" 


<690 


4xl0 10 


J033237 


2.64 


l.OxlO" 18 


(571°) 


630±50 


0.33 


8.96 


451+ 69 


60001"°°° 


6700+ 1500 


9xl0 10 



Table Notes. Ha properties from VLT-ISAAC spectra of the BLAST 250/xm sample. The seven at top have secure redshifts (the six 
which have Ha properties are used in our analysis) while the bottom two have tentative redshifts (we calculate their Ha properties, 
but exclude them from aggregate property analysis in section [3] despite being illustrated in figures). FWHM has been deconvolved with 
the instrumental resolution measured from skylines in the vici nity of Ha (~6.5A in K-band and ~4.4A in H-band), and SFR is derived 
from Ljja using the relation SFR = 7.9 X 10" 35 L£j a from Kcnnicutt (1998). The metallicity measure ments ( a ) are compute d by 
(NII/Ha) = log(S N ii/S Ha ) and (O/H) = 12+log(0/H) derived from (NII/Ha) using methods described in lMaiolino et al.1 feOOSf) . The 
characteristic uncertainty on [Nil]/ Ha is ~0.10. Stellar masses ( b ) are measured from Spitzer-IRAC photometry which brackets the 
rest-1.6/xm stellar bump (see § 2); the characteristic uncertainty in stellar mass is ~2xl0 10 Ts/Iq. 



We measure giR, the ratio of integrated IR flux to radio 
flux, as described in detail by 110 to verify the FIR/radio 
correlation in our sample. 110 finds a mean qir = 2.41±0.20 
based on the larger sample of BLAST sources with pho- 
tometric redshifts. Assuming a radio synchrotron slope of 
a=0.75, we measure a mean qi r = 2.46±0.18 w hich agrees 
with 110 and earlier findings (e.g. lDale et al.ll2007T ) that there 
is no evidence for evolution in qir with redshift. Since a has 
a significant impact on the calculation of qir, we consider 
the impact of variations in a: 110 explicitly measured a for a 
high-redshift subset of their BLAST sample and found a me- 
dian value of a = 0.4. If we use a = 0.4 instead to calculate 
qm, we measure qm = 2.27±0.17. This is still in agreement 
with 110's measurement of the FIR-radio correlation at high 
redshift within uncertainties. 

We note that iKovacs et al.l |2006) concluded that the 
local FIR/radio correlation overestimates FIR luminosity by 
factors of ~0.2-0.4dex for SMGs, which contrasts with our 
and 110's result for BLAST sources. However, the differ- 
ence is due to different FIR SED fitting procedures; when 
we refit t he 21 SMG FIR poi nts (measured at 350/im and 
850/im) in lKovacs et al.l (|2006l ) using the methods described 
in this paper (for fixed /3 = 1.5) we find qm — 2.46±0.19. 
Figure [5] shows the flux densities of the BLAST sample in 
the FIR-radio context, plotted against two t emplate SEDs 
which follow the relation: an M82 template llBressan et al.l 
120021 ) and a composite SMG SED (|Pope et alj|2008h . Some 
sources exceed the relation (often due to flux excesses to- 
wards shorter wavelengths) while others' FIR luminosities 
are over-predicted, most likely due to AGN contribution to 
radio luminosity (e.g. J033152). However, these variations 
do not appear to be weighted in either direction, indicating 
that the FIR/radio correlation is predicting FIR luminosi- 
ties accurately for the population on a whole. 




10 1 10 2 10 3 10 4 10 5 
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Figure 5. A comparison of the BLAST FIR and radio flux 
densi ties against two tem plate SEDs: an M82 spectrum (black) 
from I IBressan et al.ll20 02) and a composite SMG spectrum (blue) 
taken from lPope et al.l l|2008h . The upper panel plots all sources' 
flux densities in their rest-frame (with overplotted best-fit FIR 
SEDs (light gray) from Fig- HJ . The lower panel plots all sources 
with fluxes renormalised to their radio points. The scatter around 
the template SEDs in the lower panel (where the flux has been 
renormalised to the radio) indicates our observed variation in the 
FIR-radio correlation. While there appears to be a large scatter, 
there is no systematic offset from the FIR-radio predicted FIR 
luminosities and those measured directly in our sample. 

3 DISCUSSION 

Combining the unique BLAST FIR data with previous 
Spttzer observations and our spectroscopic redshift iden- 
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Figure 4. The FIR SEDs for high-2 BLAST sources, including data from Spitzer-MIPS (24pm, and 70pm, 160pm where available), 
BLAST (250pm, 350pm, and 500pm) and LABOCA (870pm). We fit three SEDs to the data: (1) a single dust temperature with fixed 
= 2 modified blackbody as in Eq. [T] (das hed-dotted l i ne), ( 2) a modified blackbody with treated as a free parameter (solid black 
line), and (3) an SMG composite SED from lPope et all (2008) normalised to 24pm flux density (thin gray line). The galaxies' names, 
redshifts, best fit dust temperatures, and FIR, luminosities are inset on each SED plot. The two galaxies marked TENTATIVE (from their 
tentative redshift identifications) are also the two galaxies whose fixed = 2 SED fit was significantly better than the beta-free model. 
Neither tentative galaxies are included in the analysis of section [3] 



tifications allows a full characterisation of FIR SEDs for 
250pm-bright HyLIRGs. Since 250pm mapping has the ad- 
vantage of sampling blackbody emission near its peak at 
z~2, 250pm ULIRG/HyLIRG selection is far less biased to- 
wards certain SED shapes, making 250pm-bright galaxies a 
much cleaner, unbiased subset all high-redshift ULIRGs. 



3.1 Source Density 

The seven galaxies with secure redshifts suggest that these 
250m-bright galaxies have roughly 1/5 the volume density 
of similarly luminous SMGs. We exclude the sources which 
only have tentative redshift identifications from this anal- 
ysis, although we note that they would not skew or affect 
the interpretations we draw from the entire sample. We es- 
timate this lower limit to P250 for our sample by taking the 
redshift ranges z = 1.3 — 1.7 and z = 2.0 — 2.6, given the 
gap in the atmospheric opacity between H and if -bands. 
We treat the D10 and 110 sample separately as the rela- 
tive selection depths differ (S250 >33mJy and S250 >59mJy 



respectively), finding > 8 x 10~ 6 h 3 Mpc -3 and > 5 x 
10 -6 h 3 Mpc -3 respectively. Of the spectr oscopically iden- 
tified SMG samples in the literature (e.g. IChapman et al.l 
120051 ). 51% of the sources have L FIR > 8 x 10 12 L (a 
cutoff corresponding approximately to the BLAST 250pm 
depth), which implies a luminosity-limited volume density 
of 2.5 x 10~ 5 h 3 Mpc -3 for SMGs. As much deeper, more 
uniform 250pm data become available from Herschel, the 
overlap w ith the SMG p opulation is being explored more 
fully (e.g. lElbaz et~aT]|201oh . 

Of our nine spectr oscopic sour c es, on ly four have been 
detected as SMGs in IWeifi et all (|2009h . While 250pm- 
bright sources at high-z are more rare than SMGs, the 
fact that 55% (5/9) of our sample are submm- faint (with 
(Ss7o) ~ 2mJy) highlights that the SMG population rep- 
resents only a sub s et of high-redsh i ft UL IRG activity, as 
ICasev et all (|2009T ), [chapman et all l|2004 ) and lBlain et all 
I 2004 ) suggest. The addition of 250pm-selected, submm- 
faint galaxies to the previously-known HyLIRG popula- 
tion could imply that the volume density of known high-z 
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HyLIRGs would increase from the SMG estimate roughly 
by 12%. However, more spectral observations of simi- 
lar 250^m-bright objects from Herschel are needed to boost 
these statistics and understand the actual level of contribu- 
tion. 

3.2 Near-infrared Spectral Features 

The Ha derived star formation rates of the BLAST HyLIRG 
sample underestimate the FIR SFRs by ~10x, as is of- 
ten the case with rest-UV or rest-optical emission line 
star-formation indicators in dust-obscured starburst galax- 
ies. Howe ver, we note near-IR spectroscopic observations of 
SMGs in iTakata et ail (2006) measured internal extinction 
factors of Av — 2.9±0.5 using Hq/H/3 ratios. When cor- 
recting the Ha-inferred SFRs in Table [2] for this dust ex- 
tinction the FIR-inferred SFRs are recovered, averaging to 
~2000 Me yr" 1 . This indicates that dust obscuration is sig- 
nificant in the near-IR and must be corrected for to under- 
stand the true nature of the ultraluminous activity in these 
galaxies. 

Placing our [Nil]/ Ha metallicity measurements in a 
larger galaxy evolution context, the metallicities of this 
sample (measured by converting to (O/H), i.e. (12 + 
log(0/H)} = 8.65), agree within uncertainties with the ob- 
served metallicities of the mo st massive z ^2 ga laxies, (12 + 
log(0/H)) ~8.55±0.07, in lErb et all (|2006l ). The mean 
[Nil] / Ha ratio fo r this sample, 0.29±0.23 , agrees within un- 
certainty with the lSwinbank et alj (|2004 ) SMGs, 0.41±0.38. 
While evolutionary conclusions should not be drawn from 
these data alone, the results are consistent with conjecture 
that the ULIRG phenomenon occurs at the early stages of 
a burst in star formation triggered by the merger of two 
typical gas-rich massive galaxies at z~2. 

3.3 Temperature Fitting and Selection 

Figure [6] shows dust temperature (Tdust) against FIR lumi- 
nosity, with BLAST 250/im sources and SMGs overplotted. 
Representative la detection boundaries at 70/im, 250/im 
and 870/im are shown to illustrate the populations' se- 
lection biases (2a, 3<t or 5a detection limits would have 
the same shape but be shifted to the right in luminos- 
ity; e.g. the 3a detection limit corresponds to a luminos- 
ity shift of ~0.4dex). The mean dust temperature of our 
sample, when fit with single modified blackbody SEDs, 
is 52±6K, which is comparable to the mean dust tem- 
perature of local ULIRGs of similar ( ^ 10 13 LfT] ) lum i- 
nosities, 45L10K l|Chapman et al.ll2003l ; iRieke et alj|200Sft . 
and only ~5K warmer than SMGs of similar luminosities 
>10 13 L (andis 15 K warmer than SMGs on average, which 
are 36±7K). Overall, all BLAST sources (except the lower- 
redshift J033246) have dust temperatures consistent with 
the high luminosity end of the SMG distribution. 

It is important to note that the SMG FIR luminosi- 
ties shown here are derived from radio luminosity, via the 
FIR/radio correlation, and that the associated dust temper- 
ature fits are reliant upon that assumption. To first order, 
we and others have shown that the FIR/radio correlation 
holds at these redshifts and luminosities (for direct compar- 
ison see Table [2]), however scatter is significant, with differ- 
ences in measured/derived FIR luminosities of ±0.7 dex. As 



is oft en done for literature SMGs to date l|Chapman et al.l 
2005), dust temperature is measured by using a single FIR 
data point (e.g. observed 850/im) and forcing an SED with 
fixed radio-inferred I/fir- If the radio-inferred FIR luminos- 
ity is significantly different from the actual FIR luminosity, 
then the dust temperature will either be grossly over- or 
underestimated. 

We use the BLAST sample and its full SED information 
(thus directly measuring FIR luminosity) to test the accu- 
racy of FIR fits and derived dust temperature for SMGs and 
other high-z ULIRGs. Regardless of the accuracy of radio- 
derived FIR luminosities, we find that Tdust is systematically 
underestimated by 12±19K when derived from 870/im flux 
densities. Similarly, we also measure du st temperature from 
the 70/im flux densities (as is done in ICasev et al. 2009c, 
for 70/mi- luminous radio galaxies), and find that they are 
overestimated systematically by 6±10K. 

The severity of these over- and under predicted dust 
temperatures is due in part to the sample selection. Because 
the sample is selected at 250/im, it is likely that consid- 
ering only the 870/mi or 70/im points will produce larger 
Td error than the 870/im or 70/im selected samples, simply 
because of the temperature-weighting and biasing of these 
selection wavelengths. In other words, if a galaxy is 870/im- 
bright and is 250/im-faint, then it is far more likely to have 
a cooler inherent temperature than a galaxy which is bright 
at both wavelengths. This highlights the difficulty with fit- 
ting dust temperatures to single FIR flux measurements and 
demonstrates that the luminosity-temperature distribution 
of previously studied ULIRG populations should be revis- 
ited when more complete S ED information is gat hered from 
Herschel and Scuba-2 (e.g. iMagnelli et al.|[201ol '). 

3.4 HyLIRG Evolution 

Figure [5] highlights where 250/im observations are more sen- 
sitive to hotter dust temperatures than ~870/im at z > 1. 
The sparsity of detections in the highlighted region of Fig. [6] 
indicates that hot-dust HyLIRGs are genuinely much more 
rare than their cold-dust analogues at Lfir > 10 13 L at 
high redshift. 

The dearth of hot-dust ULIRGs ( > 60 K) from these 
data is only significant in the HyLIRG (>10 13 L ) regime 
for redshifts above z = 1.5 (in other words, it is also signifi- 
cant at fainter luminosities at lower redshifts but not fainter 
luminosities at high redshifts). Due to the 250/im BLAST 
sensitivity, we have only one 10 12 < Lfir < 10 13 L ULIRG 
in our sample, and its redshift is ~1.3. This leaves the pos- 
sibility that z^2 hotter-dust U LIRGs exist, but lie beneath 
current 250/im imaging depth. ICasev et al.l (|2009ch showed 
that at slightly lower redshifts, z~1.5, star-formation dom- 
inated hot-dust ULIRGs (T d ~ 52 K, L F ir £ 2xl0 12 L ) 
have been observed at 1/5 the volume density of SMGs, but 
limitations in Spitzer-MIPS 70/im dep th prevented signifi - 
cant detections at z ^ 2. Also, work by Casey et al.l (|2009a) 
argues that hotter-dust >60K HyLIRGs at z~2 are less 
prevalent based on CO observations of submm-faint radio 
galaxies. After accounting for selection bias, the submm- 
faint ULIRG sample was ~2x less luminous in Lfir and 
Leo than CO-observed, cold-dust SMGs. 

If we assume a priori that high-z ULIRGs have the 
same dust temperature distribution as local ULIRGs (which 
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Figure 6. FIR luminosity against dust temperature for the BLAST 250/im sample (blue circles, labelled by the last three digits in 
their right ascension). The dashed (z = 2) and dotted (z = 1) lines indicate rough 1-cr boundaries at 250/^m (blue. <T250 = llmJy, 
BLAST), 70/xm (red, <T7o = lmJy, MIPS), and 870/im (black, crg7o=3mJy, LABOCA), where sources at the given redshift would have 1-cr 
significance if it has Lfir and Td corresponding to the boundary, and would be >1<t if it lies to the right of the boundary. To translate 
these curves into 3-cr detection limits, they would be shifted ~0.4 dex to the rig ht in luminosity space; the shape of the boundary curves 
woul d be maintained. 70Mm- brigh t "hot-dust" ULIRG s from lCasev et al.l j2009d ) are overlaid as red triangles. The 850^tm-detected SMGs 
from IChapman et aL I d2005h and iKovacs et al. are overlaid as small black crosses. The area enclosed by the green highlights a 

phase-space where 250/^m observations are more sensitive than 870/^m at all redshifts z > 1. 



have (T d ) =45±10K above 10 13 K), then there is a ~60% 
chance that no >60K sources are detected within a sam- 
ple of seven sources (given a Gaussian distribution of dust 
temperatures for systems >10 13 L©). This illustrates how 
limiting our sample size is when it comes to drawing con- 
clusions for the whole 250/im-luminous population. For ex- 
ample, a sample of ~30 sources with Td ust < 60 K must be 
detected in order for that likelihood to drop to £ 13%. More 
spectroscopic observations and FIR characterisations of sim- 
ilar samples are needed from Herschel and SCUBA-2 for real 
progress to be made in high- z ULIRG evolutionary studies, 
and to probe the differences with local ULIRG populations. 

Despite its significant uncertainty given the small sam- 
ple size, the lack of hot-dust systems (>60K) in the 
BLAST HyLIRG sample is consistent with predictions 
from SPH simulations for infrared luminous galaxies (see 
iNaravanan et all 2009). They suggest that the brightest 
high-z starbursts (>10 13 Lq) are at their most active phase 
during the early stages of final merger in-fall, when gas and 
dust are diffuse, extended, and cold. Warmer dust is sug- 



gested to condense either at a later stage merger, when gas 
and dust has collapsed and heated, or when it has been 
heated by a growing AGN. While this sample of galaxies ex- 
hibits warm dust (30< Td ust <50K), we have not detected 
any >10 13 Lq hot-dust (>60K) systems in this sample and 
we find a modest AGN fraction (20%); therefore, our results 
loosely support the theory that the most luminous HyLIRGs 
are triggered by major merging events. 



3.5 Confusion 

It is important to once again consider the impact of confu- 
sion limitations and deboosting factors on our conclusions. 
We excluded the uncertainty in the deboosting factor from 
our results (as discussed in § [2]) because its blind propaga- 
tion into all of the bands in our SED fits is not justified. 
As is discussed, the deboosting uncertainty is likely to have 
far greater effect on Lfir than on Td or f3 since boosting is 
correlated between FIR bands. Using a naive model where 
the main contribution to source confusion is a single addi- 
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tional source within the beam, whose FIR colours have the 
same distribution as seen in Herschel populations, we esti- 
mate that the derived Lfir is uncertain by ~0.1dex and 
Td us t is uncertain by ~9K. However, without a proper un- 
derstanding of the sources which boost the flux of our high-z 
BLAST sample, it is very difficult to determine how the dust 
temperatures might change, although it is unlikely that the 
mean would shift far from the current mean, 52±6K. This 
differential boosting issue should be investigated carefully 
with future, large Herschel 250/im-selected samples. 

The deboosting effect on luminosity is easier to quan- 
tify than the effect on Td- Th e maximum uncer tainty for 
the deboosting factor found bv lEales et all i|2009l ) is ~50%, 
which would propagate to a factor of ~2 in luminosity. The 
mean luminosity of our sample is ~3xl0 13 L©, a factor of 
~5-10 greater than most high- z ULIRG populations in the 
literature. The factor ~2 difference caused by potential de- 
boosting corrections is not significant in comparison. These 
250/im-bright galaxies are still "HyLIRGs," thus amongst 
the most luminous, extreme starbursts measured at high- 
redshift. 



4 CONCLUSIONS 

The redshift identification of these 250pim-bright, z~2 
HyLIRGs has allowed a characterisation of their near-IR and 
FIR properties, leading to the following conclusions: 

Near-IR spectroscopy (as probed here by VLT ISAAC) 
is an efficient way of identifying redshifts (50-70% success 
rates) for FIR sources which have secure photometric red- 
shifts. The redshift range of our sample is z =1.3-2.6, av- 
eraging to (z) =2.0±0.4, making up the z ^ 1 subset of 
250/im-bright BLAST galaxies. We also find that Ha star 
formation rates underpredict their FIR SFRs by ~35 times, 
and we measure metallicities which are in agreement with 
other high-2 galaxy samples, include those of much lower 
luminosities but of similar stellar mass. 

Having multiple FIR flux densities available for each ob- 
ject, we fit FIR blackbody SEDs to each source and constrain 
Lfik, Td, and f3 independent of radio flux density or mid-IR 
flux densities. We find that the FIR/radio correlation holds, 
but that SMG composite spectra, when fit to 24/xm flux 
densities, do not successfully describe the FIR properties of 
this sample. We measure FIR luminosities ~3x 10 w L©, and 
dust temperatures averaging Td us t = 52±6 K. However, we 
warn that both of these conclusions are sensitive to the ef- 
fects of flux boosting in the FIR, although we estimate that 
this should not change Lfir, by more than a factor of -~2x 
and Tdust beyond its quoted error. 

Since 250/xm selection is more sensitive to the detec- 
tion of hotter dust sources than SMG selection (at 850/im), 
the lack of >60 K hot-dust galaxies in our sample is poten- 
tially an indication that high-z, high-L galaxies are more 
extended (with diffuse, cool dust) on a whole than local 
ULIRGs. However our small sample size limits this conclu- 
sion to only < 40% likelihood. A lack of >60 K hot-dust 
specimens in the BLAST 250/im population could be telling 
to the galaxies' evolutionary stage; this work highlights the 
need for more observations of larger samples of similar and 
fainter sources. FIR mapping at 70-500/^m from Herschel 
and SCUBA2 will further select rare and poorly studied high- 



z ULIRG populations like the galaxies presented here, and 
near-IR spectroscopic observations will enable further red- 
shift identification of their counterparts, leading to a better 
characterisation of the ULIRG phenomenon at high-z. 
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